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Abstract 

POLDIS is a Monte Carlo program for polarized (semi-inclusive) deep inelastic scat- 
tering (DIS). Unpolarized DIS events are generated with the existing lepto-production 
event generators LEPTO for DIS and AROMA for Heavy Flavor production. The relevant 
spin asymmetries are computed at partonic level to first order in a s for each generated 
event, and are then convoluted with the corresponding ratio between the polarized and 
unpolarized parton distribution functions (i.e. parton polarization). This procedure pro- 
vides a polarization weight for each event. The average of these polarization weights gives 
the polarized cross section spin-asymmetry for the generated sample. The code consists of 
a set of subroutines to be linked with LEPTO and/or AROMA. No modification to these pro- 
grams is required. Some existing polarized parton distribution functions are also included. 
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1 Introduction 



Polarized lepton-nucleon deep inelastic scattering (DIS) has been studied in the last decades by 
several experiments which have measured spin asymmetries over a wide kinematic range [|l], |2j . 
These experiments have determined the spin structure functions of the proton and the neutron 
and have tested the related sum rules. When interpreted in the framework of the quark-parton 
model the experimental results show that the quark spins account for only a rather small fraction 
of the nucleon spin, thus implying an appreciable contribution either of gluons or possibly of 
orbital angular momentum. These data indicate also a large positive contribution of u quarks, 
a negative contribution of d quarks, and, surprisingly, a small negative contribution of s quarks 
to the proton spin JTJ. A general introduction to this subject can be found for example in |J. 

Inclusive polarized DIS measurements do not allow one to distinguish the role of each indi- 
vidual partonic component. A further separation of the contributions of different constituents 
to the nucleon spin, like As or AG, requires additional input from the study of semi-inclusive 
DIS, for which only limited data have been obtained so far ||. In these experiments, in addi- 
tion to the scattered lepton, one detects also one or more hadrons produced in the interaction. 
For instance, the study of polarized open-charm lepto-production allows to access the gluon 
polarization AG in a polarized nucleon ||. 

In this work we present POLDIS, a program designed to simulate polarized DIS experi- 
ments, with particular emphasis on semi-inclusive DIS and Heavy Flavor lepto-production in 
the quasi- real photo-production limit (Q 2 — > 0). For these processes, POLDIS generates the 
spin-dependent cross section asymmetries between parallel and antiparallel configurations of 
the incident lepton beam and target nucleon (or proton beam) polarizations, which are mea- 
sured in these experiments. Since these asymmetries are ratios of cross sections no absolute 
normalization is needed in their evaluation. The spin-dependent cross sections can be extracted 
from these spin asymmetries and the spin-independent cross section. 

In this program electromagnetic processes mediated by one photon exchange are imple- 
mented. The present code can be used over a wide kinematical range, where the effects of the 
weak interaction can be neglected, i.e. for Q 2 < 100 GeV 2 <C M| . 

The implementation of the polarization in POLDIS can be summarized in the following steps: 

1 - generation of an unpolarized event, 

2 - calculation of the partonic level hard-scattering spin asymmetry for this event, 

3 - evaluation of the final spin asymmetry and of the spin-dependent cross sections. 

The unpolarized event generation is performed with the LEPTO |7[] Monte Carlo and the 
AROMA || code for Heavy Flavor (HF) production. The hadronization is based on the LUND 
string model, which is known to reproduce fairly accurately the final hadronic state in a vari- 
ety of processes, and is performed with JETSET |J. The hard-scattering spin asymmetries are 
calculated for each generated event to order and are convoluted with the ratio between the 
corresponding polarized and unpolarized parton densities {i.e. parton polarization). A polar- 
ization asymmetry weight is thus obtained, and the average of these weights for the generated 
sample gives the polarized cross section asymmetry. The spin-dependent cross sections can be 
obtained from this asymmetry and the spin-independent cross section. These calculations are 
performed in a set of subroutines to be linked with the existing unpolarized lepto-production 
event generator LEPTO @ for DIS and AROMA [§ for HF production, and JETSET for the 
hadronization. No modification to these programs is required. The unpolarized parton densi- 
ties are obtained from the PDFLIB library [fLOfl ; various polarized parton densities can be selected 



1 At present, the spin dependent hard cross sections are calculated to order a s only. 



from a collection of some existing parametrizations provided with this program, although in a 
less standardized form. POLDIS has been tested with the most recent versions of these programs 
(LEPTO 6.5, AROMA 2.2, and JETSET 7.4); it is also backward compatible with older versions. 
A Monte Carlo code for polarized DIS with similar aims, PEPSI, has been presented a few 



years ago |TTfl. Similar results could, in principle, be obtained with PEPSI; however no Heavy 
Flavor generation is included in that program. In order to generate the spin asymmetries, which 
are indeed the measured quantities, it requires separate runs for opposite spin configuration in 
addition to a run without polarization. This results in a less convenient usage compared to that 
adopted in POLDIS. Additionally, POLDIS can also generate simultaneously different asymmetry 
values using different polarized parton densities. 

In the next Section we present the kinematics, formalism, and formulae for the polarized 
DIS. The partonic level hard-scattering spin- independent and spin-dependent cross sections, 
calculated to first order in a s are summarized in Appendix A and B. Section 3 describes the 
structure of the program and the implementation of the physics presented in Section 2. The 
usage of the program is explained in Section 4. We conclude in Section 5 with a comparison 
between the scattering asymmetries simulated by our program for some reactions and experi- 
mental data. 



2 Polarized cross sections 



2.1 Kinematics 

Figure [I] depicts a deep inelastic scattering (DIS) event. The four-vectors = (E, k) and 
k'^ = (E', k') represent the momenta of the incoming and scattered lepton, respectively, and 
g M = k 11 — k ,fl is the momentum transfer from the lepton to the hadron (7* four- momentum) . 
The target nucleon of mass M has four-momentum p M , and pf is the four-momentum of the 
i th hadron produced in the interaction. The interaction is usually described with the following 
variables: 

Q 2 = -q 2 « 2EE'{l-cos6) 

v = p-q/M = E-E' 

(1) 

y — p ■ q/p • k — v/E 

x = Q 2 /2p-q = Q 2 /2Mv. 

The right-hand side of each equation is valid only in the laboratory frame, where the target 
nucleon is at rest, and 9 is the lepton scattering angle. The kinematics for inclusive scattering, 
integrated in azimuth, is completely described by two of the four variables given above, say x 
and Q 2 . When describing semi-inclusive scattering, three additional variables are needed for 
each measured hadron; a common choice for these variables is the energy fraction of the hadron 
with respect to the 7* energy 

z { = Pi- p/p ■ q = Ei/is , (2) 

the hadron transverse momentum px with respect to the 7* direction, and the azimuthal angle 
between the scattered lepton and the outgoing hadron. 



2.2 Cross section asymmetries 

In a polarized DIS experiment one measures the asymmetries 



da^ — da^ da^~* — da^ 

A{1 = dan + da n and ^ = da i^ + da ^ ( 3 ) 



Figure 1: Deep inelastic scattering event. 



for longitudinal and transverse configurations of the incident lepton and target polarizations. 
The spin orientations in Eq. prefer to the laboratory frame, where the target nucleon is at rest. 
These asymmetries are directly related to the polarized structure functions g± and g 2 . In this 
paper only the longitudinal asymmetry An is discussed (and included in POLDIS). In the next 
pages, we will use the following notation for this asymmetry: An = An . 

Usually the scattering asymmetry results are presented in terms of the virtual photon asym- 
metries A\ and A 2 

A\ = — ^ ^ and A 2 = — — (4) 

°l/2 + °3/2 a l/2 + a 3/2 

where <rj is the virtual photon absorption cross section in a configuration with total angular 
momentum J along the incident photon direction, and a TL is the interference term between 
transverse and longitudinal virtual photon nucleon scattering. 
A\ and A 2 are related to the measured asymmetry All by 

A LL = D(A 1 +r ] A 2 ) (5) 

where 



^ y(y-2) ^ 2(1 -y)yjg 

~ y 2 + 2(l-y)(l + J2) ' V ~ y(2-y) E 1 ) 

in the high energy limit (large v) and neglecting the incident lepton mass. D is the depolariza- 
tion factor of the virtual photon with respect to the incident lepton, and 

R=^ (7) 



a 



is the ratio between the unpolarized cross section for the longitudinal and transverse virtual 
photon components. The ratio R = R(x, Q 2 ) can be obtained from the QCD analysis of 
unpolarized inclusive DIS data. In practice, however, one uses parametrizations of R obtained 
directly from DIS experiments (see Section p.3| ). 



The virtual photon asymmetries are bounded by the positivity relations 

\Ai\ < 1 and \A 2 (x)\ < \]R{x) . (8) 

Since also r\ -C 1 in the kinematic range of most high energy experiments, the term proportional 
to A 2 can be neglected, and 

A, - ■ (9) 
In POLDIS both asymmetries, A L l and Ai, are generated. 




a) 




b) 




Figure 2: Lowest order Feynman diagrams for DIS: a) leading order, b) Compton, c) Photon- 
gluon fusion. 



2.3 Partonic cross sections 

Owing to factorization, the unpolarized (polarized) DIS cross section can be written as a con- 
volution of the unpolarized (polarized) parton distribution function F (AF) with the partonic 
hard-scattering cross sections da (dAa): 



da A ~ F ® da + XAF <g> dAa . 

Here A refers to the parallel ft (^ = +1) an d antiparallel f I (A : 
incoming lepton and target nucleon in the 7* — N cm.. 
The terms da and dAa are the spin-independent 



(10) 

•1) spin configuration of the 



da 



- (da^ 
2 v 



da n : 



and spin-dependent 



dAa = - (da TT - da n ) 



'12) 



parts of the partonic hard cross section. One introduces also the partonic asymmetry Ull for 
the hard-scattering process 

dAa fin 

°>LL = —77- ■ (13) 

da 



The scattering asymmetry All for the reaction is obtained from 



.4 



LL 



E/dAaAF 
E/daF 



(14) 



where the sum runs over the hard-scattering sub-processes calculated to first order in a s (the 
corresponding Feynman diagrams for da are shown in Fig. ^]), and is integrated over the acces- 
sible phase space. Using the partonic scattering asymmetry Ull-, the asymmetry in Eq. can 
be rewritten, as 

E/daF a LL AF/F 
= ETdaF (15) 

In the Monte Carlo calculation, da F can be viewed as the cross section weight, and 
Ull AF / F as the asymmetry weight for the generated event. The integrals in Eq. |1^ are 
performed by summing these quantities for all events in the generated sample. 



2.4 Hard cross sections at order a s in QCD 



The leading order (L.O.) parton level process is the virtual photo-absorption 7* + q —>■ q 
(Fig. §a). At first order in QCD the gluon radiation (Compton diagram) 7* + q — > q + G 
(Fig. ^b) and the photon-gluon fusion (PGF) 7* + G — > q + q (Fig. |2]c) also contribute to the 
DIS cross section. Since for the latter two processes there are two partons in the final state, 
the first order matrix elements involve three new degrees of freedom in addition to the two 
variables, say x and Q 2 , needed for the L.O. DIS diagram (Fig. |2]a). These three new degrees of 
freedom correspond to the energy, polar angle, and the azimuthal angle between the lepton 
and QCD scattering plane of one of the final partons (the other is fixed by the kinematics). A 
suitable choice for these new degrees of freedom is |I2| : 



X P-Pq , (pXl)-(pXfq) na , 

x p — - , z Q — , <p = - — (16) 

S P-Q \pxl\\pxp q \ 

where £ is the momentum fraction of the incoming parton, and p q is the momentum of the final 
quark, and the cross sections are five-fold differential 



d 5 <r(x, Q 



2 



, Xp, Zq, ' 



dx dQ 2 dx p dz q * 



(17) 



In the virtual boson-parton cm. frame the unpolarized cross section da can be decomposed 
as0: 

da = dcr + cos (ft d&\ + cos 2<p da 2 (18) 
(note d&i = d&i(x, Q 2 , x p , z q )) and the polarized cross section dAcr as: 

dAa = dA<x + cos <p dAai . (19) 

The cos 2(p term does not appear in Eq. [19|, because it enters only in the cross section for the 
virtual photon longitudinal component, and therefore cancels in Act. After integration over the 



azimuthal angle 0, only the first term on the right-hand side of Eqs. [18] and ^ remains (i.e. 
d<j and dAcr ). 

When studying HF production, the masses of the quarks must be taken into account. In 
this case the helicity does not coincide with the quark spin: q±\/2 = qn/L + O (m/\f£\ (the 
subscript ±1/2 denotes the quark helicity, and q^/L = 1/2(1 ±75)0). The HF photon-gluon 
fusion spin asymmetry a^L reaches the value —1 of the massless case only in the asymptotic 
limit of very high energies, while at threshold Ull — +1- In the HF lepto-production via the 
PGF the relevent scales of the process are set by the HF quark mass mq, and, therefore, this 
process can be studied also in the quasi-real photo-production limit of Q 2 — » 0. 

In Appendix A we summarize these partonic cross sections calculated for one photon ex- 
change following the cross section decomposition of Eqs. [18] and pi]. In Appendix B we sum- 
marize the same cross sections expressed in terms of the Mandelstam variables s, t, and u 



integrated over the azimuthal angle <p. The unpolarized cross sections have been derived in [12 



for the massless case and in [I3| for HF. The polarized ones were re-derived by us, extending 



also the results of |TT], [L4|, |15 |. 



The partonic scattering asymmetry a LL is obtained from Eq. [13| by adding up the various 



terms of the cross sections in Eqs. 18 and 19, which are summarized in Appendix A. For 



instance, the L.O. scattering asymmetry is 



*IT - TT | r rf ■ (20) 



POLINI (S) 
POLASYM (S) 



POLSTR (S) 

POLINTL (S) 
POLEND (S) 
ALLQ (F) 
ALLQG (F) 
ALLQQ (F) 
ALLQQHF (F) 
RPAR (F) 
POLDISU (C) 



initializes the (un) polarized parton distribution functions 
calculates the polarization weight for each generated event with calls 
to the functions for the calculation of the spin asymmetries calls 
to the subroutines for the extraction of unpolarized and polarized 
parton densities, and the calculation of R 

returns the values of the polarized parton densities at given x and 

Q 2 (PDG flavor code convention) 

contains the internal set of polarized parton densities 

gives the spin-dependent cross sections (Monte Carlo estimate) 

calculates a LL for 7 + q — > q 

calculates a LL for 7 + q — > q + G 

calculates Ull for 7 + G — > q + q 

calculates Ull for 7 + G ^ Q + <3 

gives the value of R at given x and Q 2 

contains the POLDIS settings and parameters which include also the 
asymmetry values 



Table 1: Relevant POLDIS program components: subroutines (S), functions (F), and common 
blocks (C). 

For the other Feynman diagrams shown in Fig. |2|, the partonic asymmetries result in much 
more complicated expressions. In Figure |3] we plot the partonic scattering asymmetries an 
of order function of the cm. scattering angle •&* between the incoming and outgoing 

partons for various values of the 7* momentum transfer Q 2 and of the cm. energy s. To be 
noted the Q 2 dependence of these asymmetries. The angle #* is given by 

COB&* = \ - 2z q . (21) 



3 Structure of the program 

POLDIS consists of a set of subroutines for the handling of the polarization, which are linked 
with the unpolarized lepto-production event generetor LEPTO J7[ or AROMA |J. The hadroniza- 
tion is performed with JETSET || using the LUND string model. No modifications to these 
event generators is required. The PDFLIB library [|l(J is used as a source for the unpolarized 
parton distribution functions. We assume familiarity with these programs. Various polarized 
parton distribution functions, obtained from their authors, are also included, although in a less 
standardized form. 

The general structure of POLDIS is similar to a typical Monte Carlo program using LEPTO or 
AROMA with the addition of calls to some subroutines for the polarization calculations. POLDIS 
produces in output, in addition to the standard LEPTO output, the polarized scattering asymme- 
tries An and Ai, and the spin-dependent cross sections. POLDIS (as LEPTO) is a slave program, 
in the sense that the main steering code for the administration of the event generation and the 
subsequent analysis of these events, has to be provided by the user. 

The various relevant POLDIS program components are summarized in Tab. [1|. Most of names 
start with POL. These are the only components that may be accessed by the user. The program 
settings and parameters are listed in Tab. |2|. These parameters contains the values of different 
spin asymmetries, which need to be accessed by the user. 
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Figure 3: The scattering asymmetry Ull for 7* + q — > q + G, 7* + G — > q + q, and 7* + G — > Q + Q 
(HF) as a function of the cm. scattering angle $* at fixed y = 0.7 for different values of Q 2 
and for two values of the cm. energy s. 



POLLST(l) 


polarized parton distribution function 


POLLST(2,3) 


unpolarized parton distribution function (in PDFLIB format) 


POLLST(4) 


parametrization of R 


POLLST(5-9) 


unused at present 


POLLST(IO) 


number of generated events used in the asymmetry calculation 


POLPAR(l) 


Ull for current event 


POLPAR(2) 


R (o~ L /o~ T ) for current event 


POLPAR(3) 


D (depolarization) for current event 


POLPAR(4-6) 


AF/F for current event for the subsets 


POLPAR(ll-13) 


All for current event 


POLPAR(14-16) 


A\ for current event 


POLPAR(17-19) 


All for the generated sample 


POLPAR(20-22) 


A\ for the generated sample 


POLPAR(23-28) 


unused at present 


POLPAR(29-31) 


in pb - Monte Carlo estimate associated with generated event sample 


POLPAR(32-34) 


a'^ in pb 


POLPAR(35-40) 


unused at present 



Table 2: POLDIS parameters in common block /poldisu/: POLLST is an array of integers, 
and POLPAR an array of double precision real numbers. 

• At the initialization stage standard LEPTO and/or AROMA parameters and switches are 
selected. Additionally, a parametrization for the polarized distribution functions and 
for R = a L /a T are chosen, and a kinematical interval is defined for the simulation. In 
the subroutine polini the selected polarized distribution functions are read from the 
corresponding ASCII file(s). Immediately after that the unpolarized event generator is 
initialized with a call to linit (LEPTO) or arinit (AROMA). 

• In the event loop unpolarized events are generated with calls to lepto (or aroma). The 

asymmetry is calculated for each event in the subroutine polasym, and the asymmetry 
results are stored in the common block /poldisu/. At this point the user can perform 
additional analysis on the generated event: for instance, he can select a binning e.g. in x 
or y for the asymmetries All and A\, or study semi-inclusive asymmetries by requiring 
in the event a tt + with z > 0.2. 

• In the ending stage the spin-dependent cross sections are estimated in the subroutine 
polend. The relevant results are printed in the form of a table. 



3.1 Asymmetry evaluation 

The scattering asymmetry a L L is calculated for each generated event according to the under- 
lying sub-process and the kinematic variables given by the unpolarized event generator. The 
unpolarized (F) parton densities are evaluated at the given x and Q 2 . The polarized (AF) 
parton densities are also evaluated at the given x and Q 2 for 2 or 3 subsets of the selected po- 
larized parton distribution functions set (see Section |3.2| ) , and the parton polarizations, AF/ F, 
are thus obtained. The polarization weight 



AF 

W L L = CLLL X 



(22) 



is finally calulated for all subsets. 

The values of the asymmetries A LL and A\ are updated event by event: 



A LL = ^A LL + ±w LL and A 1 = ^A 1 + ^. (23) 

N is the number of events generated so far, and D is the virtual photon depolarization (Eq. |6]). 
The values of the asymmetries corresponding to different AF's from the same subset are stored 
in the common block /poldisu/. The statistical accuracy on the asymmetries depends on the 
number of generated events N and goes as 1/y/N. To study, for instance, the x and/or Q 2 
behavior of the asymmetry, the values of All and A\ can be binned as a function of x and/or 
Q 2 . 

The spin-dependent cross sections are obtained from the unpolarized cross section a and 
the scattering asymmetry All'- 

^ = ±o°(l + A LL ) (24) 
* W = l*°(l-A LL ) (25) 

3.2 Polarized parton distribution functions 

The following polarized parton distribution functions for the proton are presently included 
(in parenthesis are shown the corresponding unpolarized parton densities, which should be 
used for consistency and are automatically selected, and the x/Q 2 range of validity of the 
parametrization; the notation adopted is that of PDFLIB): 



1 GS-95 (unpolarized: DO 1.1; range: x > 10" 5 , 4 < Q 2 < 4.5 x 10 5 GeV 2 ) 

2 GS-96LO (unpolarized: GRV-94LO; range: x > 10" 5 , 1 < Q 2 < 10 6 GeV 2 ) |7| 



3 GS-96NLO (unpolarized: MRSA'; range: x > 1(T 5 , 1< Q 2 < 10 6 GeV 2 ) |7 



4 GRSV-96LO (unpolarized GRV-94LO; range: x > 10" 4 , 0.4 < Q 2 < 10 4 GeV 2 ) |8| 



5 GRSV-96NLO (unpolarized GRV-94HOMS; range: x > 10" 4 , 0.4 < Q 2 < 10 4 GeV 2 ) pf 

Typically, each set of polarized parton densities contain two or three different parametrizations, 
obtained from the same analysis. The GS polarized parton distribution functions contain three 
subsets, referred as set A, set B, and set C. The GRSV polarized parton distribution functions 
contain two subsets, referred as standard and valence scenario. In POLDIS all subsets are used 
simultaneously, thus giving an output with two or three values for the spin asymmetries All 
and A\ corresponding to the used subsets. The parton distribution functions are stored on 
a x/Q 2 grid for each parton component and stored in ASCII files. These polarized parton 
distribution functions have been obtained from the corresponding authors |16|, [17|, [L8| . 

An internal set of polarized distribution functions is also included, mainly for debugging and 
apparatus studies. Different parametrizations can be implemented by the user by simply editing 
the subroutine polintl, which contains this internal set. This internal set can be also used, for 
instance, for evaluating the effects of a large negative sea polarization As < on the scattering 
asymmetry for a particular channel. It is assumed that these polarized parton densities scale as 
the corresponding unpolarized distribution functions (no dynamical generation of the sea and 
gluon polarization is performed). For instance, the following parametrization, based on the 
SU(6) spin structure of the proton combined with a soft gluon can be used (and it is included): 

Au v = u v - fd v , Ad v = -~d V) Aq s = 0, AG = xG 



For complex targets {i.e. containing several protons and neutrons) full isospin symmetry is 
assumed between protons and neutrons: Au£ = Ad", AG P = AG", etc., and possible nuclear 
effects are neglected. 

The selection of a polarized parton distribution function set is performed by setting the 
corresponding switch to the desired value: 1 to 5 for the polarized parton densities listed 
above, for the internal set, and —1 for no polarization. At the initialization stage in subroutine 
polini, in addition to the polarized parton density, the corresponding unpolarized one is also 
selected. Our default polarized parton distribution functions set is the GS-96LO, combined 
with the unpolarized ones of GRV-94LO. 

3.3 Parametrizations of R 

In most polarized DIS experiments the virtual photon asymmetry A\ is obtained from the mea- 
sured asymmetry All (Eq. ||, using a parametrization of R = a L ja T determined from unpolar- 
ized DIS data. A similar approach is also adopted in POLDIS. The following parametrizations 
are included (in parenthesis is shown the kinematic range over which R was estimated, and 
corresponds also to the region, where the parametrization can be used safely): 

1 SLAC (range: x > 0.03, Q 2 > 0.35 GeV 2 ) || 

2 NMC-97 (range: x > 0.003, Q 2 > 0.30 GeV 2 ) 

3 BKS-97 (range: 4 x 10~ 5 < x < 0.1, 0.01 < Q 2 < 360 GeV 2 ) 

The selection of the desired parametrization of R is obtained by setting the corresponding 
switch to the appropriate value: 1 to 3 for the parametrizations listed above, and for R = 0. 
Our default is the NMC-97 parametrization. 

3.4 POLDIS parameters 

The common block /poldisu/ contains the program settings and parameters as illustrtated in 
Tab. g 

COMMON / POLDISU / POLLST(IO), POLPAR(40) 

In this common block different asymmetry values are stored. For each spin asymmetry and 
spin-dependent cross section there are three different values, which correspond to the three 
subsets of the selected polarized parton distribution functions set. 

3.5 Interface with LEPT0 

The relevant kinematics, program parameters, settings, and switches, used also for the asym- 
metry calculations, are stored in the LEPT0 common block /leptou/: 

COMMON / LEPTOU / CUT(14), LST(40), PARL(30), X, Y, W2, Q2, U 

The correspondence between the kinematical variables stored in the common block /leptou/ 
to the ones discussed in the previous pages is: 



X = x, Q2 = Q 2 , Y = y, U = u, and 
PARL(28) = x p , PARL(29) = z q , PARL(30) = 0. 
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Figure 4: a) Simulated inclusive asymmetry A[ compared to SMC data from polarized pro- 
tons [^2| using the GS-96LO set A (full line) and the GRVS-96LO standard seen, (dashed 
line) polarized parton densities and the NMC-97 parametrization of R. b) Same as (a) for Af 



compared to SMC data from a polarized deuteron target |23 . 



Additionally, the following LEPT0 switches are used for the asymmetry evaluation: 
LST(22) specifies the struck nucleon: 1 = proton, 2 = neutron. 

LST(24) specifies the hard-scattering sub-process: 1 = q, 2 = qG, 3 = qq, 5 = QQ— event (HF). 
LST(25) specifies the struck quark: 1 = d, 2 = u, 3 = s, — 1 = d, — 2 = u, — 3 = s. 

The unpolarized cross sections measured in pb are stored in PARL(23) (numerical integra- 
tion at the initialization stage) and in PARL(24) (Monte Carlo estimate). 

The kinematics of the interaction and of all produced particles is stored in the JETSET 
common block /lujets/. 



4 How to run P0LDIS 

In addition to the standard LEPT0 (and AROMA) input parameters and switches, such as the beam 
energy, target material, etc. (we assume the user to be familiar with them), two additional 
input switches are required (see previous Section): 
POLLST(l) = to 5 for the polarized parton densities, and 
POLLST(4) = to 3 for the parametrization of R. 

As already mentioned above, the user must provide a steering code for the administration 
of the event generation and analysis. Before the initialization of P0LDIS and LEPT0 the relevant 
parameters, switches, etc. must be set to the corresponding values. 

At the end of the event generation loop POLPAR contains various asymmetry values and 
the spin-dependent cross sections associated with the subsets of the selected polarized parton 
density set. 



5 Results of test runs 

The inclusive asymmetry A\, obtained with the GS-96LO set A and the GRVS-96LO standard 
seen, polarized parton densities and the NMC-97 parametrization of R, is compared in Fig. [| 
to the SMC data from polarized proton [p2[] and deuteron p3f targets. A fairly accurate 



agreement between the simulated and the real data can be observed in these plots. It has 
to be noted, however, that polarized DIS data were used for the evaluation of the polarized 
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Figure 5: Simulated semi-inclusive asymmetry A\ + on polarized protons (a) and Af + on polar- 
ized deuterons (b) for positive hadrons with z > 0.2, compared to SMC experimental data || 
(GS-96LO set A full line, and GRVS-96LO standard seen, dashed line). 



parton distribution functions used here, and therefore this agreement can not be considered as 
a meaningful physics result. On the other hand, such an agreement shows the validity of the 
procedure adopted in the simulation and the correctness of the calculations. 

Figures [| and |6] show the semi-inclusive asymmetries for positive and negative hadrons 
generated with P0LDIS, respectively, compared to the SMC data || from polarized protons and 
deuterons. Also for this simulation we used the GS-96LO set A and the GRVS-96LO standard 
seen, polarized parton densities and the NMC-97 parametrization of R. In the Monte Carlo 
simulation a separation between charged pions and kaons can also be made. 
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Appendix A: spin-averaged and spin-dependent hard cross 
sections 

In this Appendix we summarize the various terms appearing in the unpolarized and polarized 
partonic level cross sections (Eqs. |18] and [T9| ) for the processes included in this program. To 
shorten the notation, the common factor 



2a 2 e 2 q 
~x~Q^ 



(a.i; 



describing the electromagnetic coupling is singled out. e q is the charge of the struck quark. 



7 9^9: 

For the L.O. diagram (Fig. 0a), the unpolarized cross section is: 



and the polarized one is: 



d^o = i(l + (l-y) 2 ), 



dAa = _(l-(l- y ) 2 ). 



(A.2) 



(A.3) 



Note that there is no dependence on the azimuthal angle <fi. 



For the Compton diagram (Fig. 0b) the unpolarized cross section is 



d<r 
da 2 



2a s \l + {l-yf 



3tt 
AOs 
3tt 
4a a 
3tt 



X p Z q 



+ 2 (XpZ q + 1 



(l-x p )(l-z q ) 
(V - 2) yfTTy ^—I^IZ { i -x p -z q + 2x p z q 

(1 - y) x p Zg , 



+ (1 -y) ^x p z q 



and the polarized one is: 
dAo"o = 
dA&i -- 



2a s l-{l-yf 



3tt 
3vr 



X p + Z q 



+ 2 (X p + Z q ) 



v\/i-y. 



{I - X p ){\ - Z q ) 

(1 X p Zq) . 



XpZq 



(1-X p )(l-Zq) 



(A.4) 
(A.5) 
(A.6) 

(A.7) 
(A.8) 



• 7* 9 "~ > 1 + Q (light quarks): 

For the PGF diagram (Fig. 0c) - massless case (m q = 0, q = u, d, s) the unpolarized cross 
section is: 



d<r 

d&i 
da 2 



1 + (1 - y) 2 [4 + (1 - x p ) 2 \ [z 2 q + (l-z, 



47T 



Xp(l Xp) 



\ z q {l-z q ) 



(1 - 2x p ) (1 - 2z g ) 



Oi, 



7T 



(1 -y) rr p (l-Xp) 



+ (1 - j/) 8z p (l - x p ) j> (A.9) 

(A.10) 
(A.ll) 



and the polarized one is: 
dAoo 



a s l-{l-y) 2 (2x p -l) 


^ + (1 " z q f 


4vr 2 z 9 





a. 



Xp{\ Xp) 
\ Z li l - Z q) 



(1 - 2z q ) . 



(A.12) 
(A.13) 



7* g — > Q + Q (heavy quarks): 
For the PGF diagram - massive case (mq 7^ 0, Q = c, b) the unpolarized cross section is: 



da n 



a. J 1 + (1 - 2/) 2 

47T 



xj + (1 - x p ) 2 ] + (1 - z q ) 2 ]+/3(l-2x p ) (3(2 Xp - (3) 
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(A.14) 
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(A.15) 
(A.16) 



and the polarized one is: 



a s l-(l-y) 2 

= 2 



+ 



(3 



z q (l-z q ) 2z 2 (l-z q ) 2 \ 



z 2 + (l-z q 
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d Aai = y yj 1 - y 
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(A.17) 
(A. 18) 



The factor /3 is related to the HF quark velocity v in the photon-gluon cm. by v — a/1 — JJ, 
and it is given by 



R = Ax — - 



(A.19) 



where mg is the mass of the heavy quark. 



Appendix B: spin-averaged and spin-dependent hard cross 
sections using Mandelstam variables 

In this Appendix we rewrite the unpolarized and polarized partonic cross sections (Eqs. [18] 
and |19D integrated over the azimuthal angle in terms of the Mandelstam variables s, t, and 
u. The relations between these variables and the x p , z q variables used in Appendix A are: 

§ = Qil—^L = 1Mvt-Q 2 

Xp 

u = m 2 Q - Q 2 — (B.l) 

Xp 

t = m -Q 



Xp 



§ + t + u = 2m 2 Q - Q 2 

where £ is the momentum fraction of the incident parton (as in Eq. |TB|) and rriQ is the mass of 
the (HF) quark. 

These cross sections expressed in terms of the Mandelstam variables depend on Q 2 , x, s, 
and t (d&i = d&i(x, Q 2 , s, t). The common factor describing the electromagnetic coupling 

Aiia 2 e 2 

(R2) 

is singled out as in Appendix A. The additional factor 2ir compared to Eq. |A.1| comes from the 
integration in 0. 

• 7* q -> q: 

For the L.O. diagram (Fig. 0a), the unpolarized cross section is: 

da = i(l + (l-2/) 2 ), (B.3) 

and the polarized one is: 

dAa = \{l-{l-y) 2 ). (B.4) 



l*q -> q + g: 

For the Compton diagram (Fig. ||b) the unpolarized cross section is: 
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and the polarized one is: 
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(B.5) 



(B.6) 



• 7* 9 ~ ¥ Q + Q (light quarks): 



For the PGF diagram (Fig. 0c) - massless case (m q = 0, q = u, d, s) the unpolarized cross 
section is: 



as 1 j l + (l-y) 2 Q 4 + s 2 u 2 + t 2 8Q 2 s 
An(s + Q 2 ) 2 \ 2 (s + Q 2 ) 2 ut 1 V ' (s + Q 2 ] 



and the polarized one is: 



i a - a. 1 l-{l-y) 2 Q 2 -su 2 +t 2 

dAa ° = 4^(JTW 2 JT^IT' (R8) 



1* 9 ^ Q + Q (heavy quarks): 
For the PGF diagram - massive case (mg ^ 0, Q = c, b) the unpolarized cross section is: 
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and the polarized one is: 
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